Incoherent transport induced by a single static impurity in a Heisenberg chain 
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The effect of a single static impurity on the many-body states and on the spin and thermal transport in the one- 
dimensional anisotropic Heisenberg chain at finite temperatures is studied. Whereas the pure Heisenberg model 
reveals Poisson level statistics and dissipationless transport due to integrability, we show using the numerical 
approach that a single impurity induces Wigner-Dyson level statistics and at high enough temperature incoherent 
transport within the chain, whereby the relaxation time and d.c. conductivity scale linearly with length. 



PACS numbers: 71.27.+a, 71.10.Pm, 72.10.-d 

Transport in the one-dimensional (ID) quantum system 
of interacting particles still offers several fundamental the- 
oretical challenges. In the context of the electron transport 
through barriers or weak links the role of repulsive electron- 
electron interactions within the wire was shown to be cru- 
cial |H 01 since at low temperature T — > the interaction 
within Luttinger-liquid (LL) phenomenology renormalizes the 
transmission through the barrier to zero effectively cutting the 
chain for transport yfl. Neglecting the effect of Umklapp 
processes within the wire above the Kane-Fisher temperature 
T* the transmission through the barrier should become finite 
indirectly confirmed in a numerical study of the ID spin 
model fl]. On the other hand, in the last decade it has become 
increasingly evident that the LL low-energy description is not 
enough to establish transport properties, if they are dominated 
by Umklapp processes. It has been shown ||5l 0] that pure 
ID integrable models of interacting fermions exhibit in spite 
of Umklapp at any T > dissipationless (ballistic) transport 
manifested, e.g., in a finite charge stiffness D(T > 0) > 0. 
On contrary, a generic system of interacting fermions would 
exhibit dissipation in the wire and also finite d.c. transport co- 
efficients, e.g., the d.c. conductivity a(u> — > 0) = era < oo. 
The distinction is closely linked to the statistics of many-body 
levels [7], which follow the Poisson level distribution for the 
integrable system and the Wigner-Dyson (WD) distribution 
for the generic nonintegrable system Random disorder, 
strong enough to overcome finite-size effects, in such models 
generally leads to the WD statistics for nearest levels iHH], 
while the d.c. transport seems to be normal (dissipative) at 
T>0[10], in contrast to the Anderson-type localization per- 
sisting atT = ifTHl. 

Our goal is to understand within this context the effect of a 
single static impurity in the ID integrable system of interact- 
ing particles. While of fundamental importance this question 
is also directly relevant in connection with ongoing experi- 
ments on novel quasi- ID materials where electronic proper- 
ties of the pure system can be well described within the inte- 
grable spin-1/2 Heisenberg model [12] with dilute impurities 



introduced in a controlled manner and their influence studied, 
e.g., on the thermal transport. 

In the following we show on the example of the ID spin 
model with periodic boundary conditions that a single static 
impurity can qualitatively change the level statistics to the WD 
one, noticed also in the recent study of the onset of quantum 
chaos lfl3tl although the perturbation scales as 1/L where L is 
the length of the system. At the same time, the impurity leads 
to the vanishing of the spin stiffness D at elevated T > 0. 
In such a situation, it is meaningful to discuss the decay of 
the spin and energy current within the ring and related trans- 
port rates 1/r which we show to be well defined and scale 
as 1/r oc 1/L as expected for the homogeneous wire with a 
single localized perturbed region. 

We consider the ID anisotropic Heisenberg model (AHM) 
with a single-site static impurity field, 

H = E J (Sf +1 Sf + Sf +1 Sf + AS? +1 S?) + b S%, (1) 
i 

where S a , a — x, y, z are spin-1/2 operators, J is the mag- 
netic exchange coupling (we use units J = 1), A the 
anisotropy and bo the local impurity field. Numerically we 
study chain (ring) of length L with periodic boundary condi- 
tions. 

In the absence of the impurity the AHM, Eq. £[), is inte- 
grable and as the consequence reveals the Poisson level dis- 
tribution Pp(s) = exp(— s) where s — (E n +i — E n )/Ao 
(Aq is the average level spacing) as well as the dissipation- 
less transport [g]. Let us first consider the effect of finite 
bo on the level statistics. We investigate this question by 
performing the (full) exact-diagonalization (ED) study of fi- 
nite size systems with L = 10 — 16. It should be pointed 
that in this range the number of many-body states varies (in 
the Sl ot = sector) in a wide range N st = 10 2 — 10 4 
and the corresponding Ao = 2.1CP 2 — 5.1CP 4 . The gen- 
eral conclusion is that finite &o > induces WD distribu- 
tion Pwd(s) = (7rs/2) exp(— 7rs 2 /4) following the random- 
matrix theory (RMT) |7] in spite of the fact that the perturba- 
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tion is only an 1 /L effect (the perturbation bo/ 2 relative to the 
full energy span A# ~ L J). 

To be concrete we present here two standard tests for the 
closeness of the RMT. The first one is parameter 77 |8J mea- 
suring the normalized distance to the WD distribution, 

/•so rso 

V= [P(s)-P WD {s)]ds/ [P P (s) - P WD (s)]ds, 
Jo Jo 

(2) 

where P(s) is the actual level distribution and So = 0.473 is 
chosen to be the intersection of Pp (s) and Pyy r> (s) @] . In or- 
der to stay within the regime of homogeneous density of states 
we analyze only one half of intermediate many-body states, 
as relevant for the high-T properties discussed here. ED re- 
sults for resulting 77 as a function of 60 for chosen intermediate 
A = 0.8 are presented for different L. To avoid the effect of 
higher degeneracy of levels at S^ ot = 0, 60 = presented re- 
sults in Fig. 1 are for Sf ot = 1. In the absence of impurity 
(bo = 0) we obtain r\ = 1 since P(s) = Pp(s) due to the in- 
tegrability of the pure AHM. The most important conclusion 
is that rather weak impurity bo ~ 0.2 in largest L = 16 causes 
a fast drop to 77 ~ 0, i.e., to P(s) ~ Pwd(s), whereby the 
threshold value of bo is decreasing with L so that for largest 
L = 16 reachable with ED we get P(s) ~ Pwd(s) in the 
range 0.2 < 60 < 1-5. On the other side, it is quite remark- 
able that r] starts to recover towards 77 ~ 1 again for large 
&o > 1. This can be easily explained by noting that large 
1 60 1 ^ 1 effectively cut the ring and lead to the AHM with 
open ends which is again an integrable model. 
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Figure 1 : Parameter 77 for the deviation from the WD level distribu- 
tion vs. impurity field bo for A = 0.8 and various L. 

Even stronger probe of the level statistics is the correlation 
A3 measuring the level fluctuations beyond the nearest neigh- 
bor levels 

1 r N 

A 3 = ^min^.s / [M{E) - AE - B] 2 dE, (3) 

<^ v J-N 

where J\f(E) is the integrated density of states with E — 
E/Ao [18]. A 3 should behave as A 3 ~ N/15 for Poisson 
distribution, and asymptotically as A3 ~ (lnA)/-7r 2 within 
the RMT 1211 . In Fig. 2 we present results for A3 (A) for fixed 



A = 0.8, b = 0.8 as obtained for different L = 12 — 16. 
A comparison with the result expected from the RMT shows 
that A3 (A) approaches the latter very accurately in an inter- 
val N < N*(L) with A* strongly (exponentially) increasing 
with L, while the deviation into a Poisson-like linear depen- 
dence A3 oc N appears for N > N*(L). Such a generic 
crossover has been observed also in other systems 01511 and 
one can discuss the relevance of the related crossover en- 
ergy scale e = A*Ao. Fast increase of N*(L) one can un- 
derstand by noting that the impurity perturbation being L- 
independent mixes up many-body levels [8] within the interval 
e whereby separation between many-body levels decreases as 
Ao oc exp(— L). We can estimate e ~ &o/(4£) within the 
XY (A = 0) model which gives right order of magnitude for 
observed A* in Fig. 2. More detailed analysis in analogy to 
other systems [8J is difficult due to the complicated nature of 
states at intermediate A. 
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Figure 2: Level-fluctuation parameter Aa(iV) for fixed A = 
0.8, 60 = 0.8 and different system length L. For comparison the 
RMT result is presented (dotted line). 

Closely related to the onset of the WD distribution by a 
single impurity is the vanishing of the T > coherent (bal- 
listic) transport characteristic for integrable systems J5I 01. 
The measure of the coherent component is for the spin trans- 
port the spin stiffness D(T) (equivalent to the charge stiff- 
ness for the related fermionic model). It can be defined via 
the gauge phase <fi into the spin-flip terms in Eq. ((T), as 

0)Sf +1 Sf. At finite T > the 
spin stiffness can be expressed as 

n n T 

where p n = exp(-f3e n )/Z with Z = J] n exp(-/3e„), and 
the last relation becomes an equality provided that the sus- 
ceptibility for persistent current vanishes (for finite systems 
at large enough T). On the other hand, D still depends on 
the value <f> where derivatives in Eq. (0]) are taken. For the 
XY model at A = corresponding via the Wigner-Jordan 
transformation to tight-binding noninteracting fermions with 
t = J/2 and a potential impurity eq — bo/ '2 one can establish 
the relation with the transmission through the barrier, as used 
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also in connection with the evaluation of the ID conductance 
lfl7ll at T — 0. For general T > one gets in the case of NI 
fermions and L — ► oo, 



5 = |EA(i-A)W 2 ft- 



where vj. = It sin k, fk is the Fermi function with \i 
half-filling (S z tot = 0) and 



(5) 



Ofor 



9k 



|i fe | 2 sin 2 (I^) 



1 



\t k | 2 cos 2 {L<t>y 



At 2 sin 2 k 



At 2 i 



■k + 4 



(6) 



Numerically we recover the behavior D(<f>) as follows from 
Eqs. ( 1516b for arbitrary &o as far as A — > 0. For A > the 
dependence on remains qualitatively similar, although ir- 
regular due to strong dependence on L. In the following we 
calculate D(L) for fixed cf> — ir/(2L). Results for A > are 
nontrivial for any T. Since results of full ED are best at high 
T, we restrict ourselves here to the limit j3 — > 0. It has been 
shown for the pure model that D/(3 remains finite and nontriv- 
ial in the thermodynamic limit L — > oo due to integrability of 
the model [6]. 

In Fig. 3 we show results for D//3 vs. 1/L for chosen A = 
0.8 and for four cases bo — 0, 0.5, 1, 2. It is evident from 
Fig. 3 that bo > cases are qualitatively different from the 
&o = where D scales linearly in 1/L towards a finite D//3 ~ 
0.035. On the other hand, bo > induces an exponential-like 
decay of D — > 0, at least for large enough L > L* and not 
too weak bo- This is closely related to the onset of the WD 
distribution and the effective breaking of the integrability. 
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Figure 3: High-T spin-stiffness D//3 vs. 1/L for fixed A = 0.8 and 
different &o- 

Rapid (exponential) vanishing of D(L — > oo) at T > is 
the indication that the transport is not ballistic and becomes 
incoherent (resistive) beyond the characteristic L*. In order 
to test this directly we evaluate dynamical spin conductivity 
<j(oj) as well as the related thermal conductivity k(lu), defined 



o-(cj) 



uj+L 



k(lu) = 



uj+L 



where j and je are spin and energy current, respectively, with 
corresponding susceptibilities 



(8) 



and analogous definition of Xj E jE Note that for 'normal' 
transport one expects \% ~ Xjji 1 ^ —> 0). In a nondissipative 
case, however, x% ~ Xjji 1 ^ ~ > 0)) — ZLD > 0. For fur- 
ther discussion it is convenient to introduce and analyze also 
corresponding memory functions M(oj) and N(u>), defined 
respectively as 





cr(w) = — 



A.; 



L uj + M(uj)' 



k(uj) 



iP Xj 



N(uj)' 



(9) 



The advantage of studying k(oj) is that jg is a conserved 
quantity in the pure AHM ||6t], hence N (oj) = and con- 
sequently N(oj) 7^ appears only due to bo ^ 0. On 
the other hand, j is not conserved and M(oj) — Mq{uS) is 
nontrivial even in the absence of impurities. Nevertheless, 
M(J(u> = 0) = at any T as required to obtain D{T) > 0. 
In the following we evaluate a(u>) and k(u) at T > using 
the microcanonical Lanczos method (MCLM) lfl4ll to calcu- 
late the dynamical susceptibilities, Eqs. ( 1817b . in systems with 
L = 16 — 24. Typically, Nl ~ 2000 Lanczos steps are used 
to obtain spectra with high u> resolution, so that the additional 
broadening is only 6 = 0.01. In the following we present 
merely results at ft — ► 0, which do not qualitatively change 
with T down to T ~ J. At least for larger L smooth <j{uj) 
and k(oj) then allow the evaluation of M{ui) and N(uj) via 
Eq. ©. 
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Figure 4: a) High-T results for scaled thermal conductivity k' /L vs. 
Luj for A = 0.5 and bo = 0.5 as calculated for different L = 16 — 
24. b) Extracted scaled memory function LN"(uS). Perturbation- 
theory result is also plotted (dashed curve), c) N"{uj)/b^ for A = 
0.5 and different impurity fields bo = 0.3 — 0.5, obtained at fixed 
L — 24. 

Let us start with the analysis of (real part) k'(u>), which re- 
veals a Lorentzian (Drude) form for A = 0.5, bo = 0.5 as 
shown in Fig. 4a. Moreover, results show universal size scal- 
ing as k'(Luj)/L. Hence, corresponding N"(w) also scales as 
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1/L, so that the size-independent quantity is N(u>) = LN(w) 
being quite structureless for u> < 1 as shown in Fig. 4b. It is 
also evident from Fig. 4c that N"(uj) oc 6q at least for weaker 
&o < 0.5. On the other hand, at larger b Q > 0.5 N" (oj) obtains 
a characteristic peak at low oj — > which strongly reduces the 
d.c. value k(uj = 0). I.e., on entering the regime bo > 1 the 
impurity starts to cut the ring for the d.c. transport. 

The regular behavior of N"(ui) for weaker bo < 1 gives 
support to the attempt to evaluate the memory function within 
the perturbation approach Hill using the force-force correla- 
tions, 



1 



(Xff(<* 



x° ff ) 



f=i[H,j E ], (10) 



where is the force-force dynamical susceptibility. Re- 

sults for Np(u>)/L evaluated using eigen-states obtained by 
the ED of the system L = 14 are for comparison also pre- 
sented in Fig. 4b. The correspondence is quite satisfactory. 
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Figure 5: a) Scaled spin conductivity cr'/L vs. Luj for A = 
0.5, bo = 0.6 for different L — 16 — 24, and b) corresponding 

M"(cj). 

Finally, we present in Fig. 5 also analogous scaled cr'/L vs. 
Loj and corresponding M"(uj) for fixed A = 0.5, bo = 0.6 
and different L = 16 — 24, as obtained via the ED (L = 16) 
and MCLM (L = 20, 24) at j3 -> 0. Since M (u) is nontriv- 
ial even for 6o = one can discuss possible decomposition 
M(w) = M (gj) + M(lo)/L. Results confirm that at low 
u> < 0.2 (j(uj) reveals a Lorentzian with M"(uo — > 0) scaling 
as 1/L, the only contribution in this regime coming from the 
impurity. Moreover, we notice that for fixed A < 1 and bo 
M"(0) and 7Y"(0) are quite similar in value, e.g., by compar- 
ing Fig. 4c and Fig. 5. This indicates that we are at A < 1 
close to the validity of the Wiedemann-Franz law requiring an 
unique transport relaxation rate. 

It should be pointed out that obtained incoherent trans- 
port is characterized with the relaxation times r and d.c. 
conductivities scaling linearly with L, as expected for ID 
systems with a single perturbed region. Hence, the length 
independent quantities are er(0)/L — x°jj I [LM" (lo — 
0)) = x*jj/M"(uj = 0) and the corresponding thermal one 
K (0)/L = Pxl jE /N"(u; = 0). 



In conclusion, we have shown that the transport in the con- 
sidered anisotropic Heisenberg model on the ring with a single 
static impurity is quite unique. Since both the spin (for A < 1) 
and thermal conductivity at any T > are dissipationless in 
the pure system, one can study directly the nontrivial effect 
of a single impurity on the level statistics and transport in the 
many -body quantum system. We have shown that single static 
impurity induces an incoherent transport with a well defined 
current relaxation time which scales as r oc L. This should be 
contrasted with the case of noninteracting fermions in Eq. ([5]l 
where a single impurity only reduces the stiffness D but does 
not lead to the current relaxation within the ring at any T. 
The fundamental difference seems to come from the Umk- 
lapp processes which are revived by the impurity and lead to 
the decoherence between successive scattering events on the 
impurity. In this sense it is also plausible that for a finite but 
low concentration Ci of static impurities in a chain (as relevant 
for experiments [12]) we expect that our results can be sim- 
ply generalized as 1/r oc 1/L — > Cj, as evident also from the 
lowest-order perturbation theory, Eq. ( fTOb . 

Although we studied here the AHM model, results and con- 
clusions could be plausibly generalized to ID integrable chain 
systems with periodic boundary conditions and a localized 
perturbed region. In this paper we presented only results in 
the high-T regime, still the phenomenon is expected to per- 
sist as far as the Umklapp processes are effective, i.e., for T 
above some characteristic Umklapp temperature Tjj. We can 
speculate that the LL phenomenology |lll|2|] can become effec- 
tive only for T <C Tjj, whereby there is another Kane-Fisher 
scale T* which divides the regimes of cut chain for T < T* 
and renormalized coherent transmission for T > T*. Fur- 
ther studies are needed to establish these scales for relevant 
impurities and models. In any case, such phenomena are ex- 
pected to be relevant in connection with recent experiments 
on thermal conductivity in spin chains where dilute impurities 
are introduced IU2I1 . 
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